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Systems with an 0(n) symmetrical Hamiltonian are considered in a d-dimensional slab geometry of 
macroscopic lateral extension and finite thickness L that undergo a continuous bulk phase transition 
in the limit L — ► oo. The effective forces induced by thermal fluctuations at and above the bulk 
critical temperature T Cj00 (thermodynamic Casimir effect) are investigated below the upper critical 
dimension d* = 4 by means of field-theoretic renormalization group methods for the case of periodic 
and special-special boundary conditions, where the latter correspond to the critical enhancement of 
the surface interactions on both boundary planes. As shown previously [Europhys. Lett. 75, 241 
(2006)], the zero modes that are present in Landau theory at T Ct00 make conventional RG-improved 
perturbation theory in 4 — e dimensions ill-defined. The revised expansion introduced there is utilized 
to compute the scaling functions of the excess free energy and the Casimir force for temperatures 
T > T c ,oo as functions of L = L/£oa, where £oo is the bulk correlation length. Scaling functions of the 
//-dependent residual free energy per area are obtained whose L — > limits are in conformity with 
previous results for the Casimir amplitudes Ac to 0(e 3 ^ 2 ) and display a more reasonable small- L 
behavior inasmuch as they approach the critical value Ac monotonically as L — > 0. Extrapolations 
to d = 3 for the Ising case n = 1 with periodic boundary conditions are in fair agreement with 
Monte Carlo results. In the case of special-special boundary conditions, extrapolations to d = 3 
are hampered by the fact that the one-loop result for the inverse finite-size susceptibility becomes 
negative for some values of L when e > 0.83. 
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I. INTRODUCTION 

When a classical or quantum fluid, or an n- vector mag- 
net with n = 1, 2, 3 is confined by macroscopic bodies 
such as two parallel plates, walls, surfaces, or interfaces, 
of area A, its free energy F depends upon the distance 
L between these boundary planes. The L dependence 
implies a force 



T C {T, L) = -k B T 



dL 



(1.1) 



between the plates, where / cx = (F/fc^T) — L/b is the 
reduced excess free energy per unit area, fb is the re- 
duced bulk free energy density, and the limit A — > oo has 
been taken [1-5]. By analogy with the familiar Casimir 
force [6] produced by vacuum fluctuations of the elec- 
tromagnetic field between two metallic plates (and slight 
abuse of language) , Tc is conventionally called "thermo- 
dynamic Casimir force." 

Provided long-range interactions are either absent or 
negligible, this force decays exponentially for separations 
L > £oo(T), where £, X (T) is the bulk correlation length. 
Near a continuous phase transition in d bulk dimensions, 
£oo(T) diverges as |T — T Ci00 | _ly at the bulk critical tem- 
perature T CiOC . Therefore, the Casimir force Tc{T,L) 
extends to distances L much larger than the microscopic 
scale a (~ radius of atoms, lattice constant). 

Writing 



f ex (T,L) = f s (T) + f ICS (T,L) , 



(1.2) 



let us decompose the excess free energy / ex into an L in- 
dependent surface part f s (T) = f cx (T, oo) and a residual 
finite-size contribution f Tes (T,L). The latter behaves as 



/res (-^c.oo ? ^) 



(1.3) 



at the bulk critical point T = T c 
the long-ranged effective force 



and hence produces 



(d-l)AcL 



(1.4) 



Here Ac , the so-called Casimir amplitude [2] , is a univer- 
sal quantity, which depends on the bulk universality class 
of the phase transition considered and gross properties of 
the boundary plates, but is independent of microscopic 
details. 

Such thermodynamic Casimir forces have been the sub- 
ject of much interest recently [7-11]. Clear experimental 
evidence for their existence has been found in the thin- 
ning of wetting layers of liquid 4 He as a function of tem- 
perature on approaching the lambda line [10, 11]. 

Near the bulk critical point the residual free energy 
density and Casimir force are expected to have the scaling 
forms 



and 



— — L t, 



(1.5) 



(1.6) 
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where and 



S(L) = (d-l)e(L)-L9'(L) 



(1.7) 



should be universal functions of L = £/£oc- These expec- 
tations rest on the assumption that £oo and L are large 
compared to other lengths, which means, in particular, 
that the symmetry breaking field h vanishes and long- 
range interactions are either absent or negligible. 

In studies of the Casimir effect in QED, matter usually 
is taken into account only through the choice of appropri- 
ate boundary conditions on the surfaces of macroscopic 
bodies. Hence they involve free field theories under given 
boundary conditions. Systematic theoretical investiga- 
tions of the Casimir effect at critical points are a much 
greater challenge in that one has to deal with interacting 
field theories in finite and bounded systems [12, 13]. 

A first fairly detailed study of the thermodynamic 
Casimir effect was made about 15 years ago by Krech and 
Dietrich (KD) for the <f> A theory on a slab R^ 1 x [0, L] 
of thickness L [7, 8]. Building on Symanzik's work [14] 
in the 80s and the simultaneously emerging field-theory 
approach to critical behavior of systems with boundaries 
[12, 13, 15-18], these authors considered five different 
boundary conditions p, namely, periodic (p = per), an- 
tiperiodic (p = ap), and the three noncquivalent combi- 
nations (D,D), (D,sp), and (sp,sp) of Dirichlet (D) and 
special (sp) boundary conditions on the slab's two bound- 
ary planes. Here the former (D) means = as usual, 
while the latter (sp) is the case of a Robin boundary con- 
dition d n 4> = ccj) for which c takes the special value c sp 
corresponding to the critical enhancement of the surface 
interactions on the respective boundary plane. 

Restricting themselves to temperatures T > T CtOQ , KD 
performed two-loop calculations for 4 — e dimensional 
slabs under these boundary conditions p and determined 
the e expansions of the Casimir amplitudes as well 
as those of the corresponding scaling functions 0^ to 
first order in e. 

In a recent paper with Shpot [19], we have shown that 
conventional renormalization-group (RG) improved per- 
turbation theory, on which both Symanzik's [14] and 
KD's [7, 8] analyses are based, becomes ill-defined at 
T c .oo beyond two-loop order due to infrared singularities 
for those boundary conditions that involve a zero mode at 
T Ci00 in Landau theory. This applies to both p = per and 
p = (sp, sp). To remedy these deficiencies, we performed 
a reorganization of field theory such that the resulting 
RG-improved perturbation theory remained meaningful 
at T c oo . It was found that the small-e expansions of the 
corresponding Casimir amplitudes A^ involve fractional 
powers t k / 2 , with k > 3, and powers of In e. Furthermore, 
explicit results for these scries to order e 3 / 2 were given. 

In this paper we will utilize this approach to compute 
the scaling functions @^ p \ and hence E^ p \ for p = per 
and (sp, sp) to the same order of RG-improved perturba- 
tion theory. The results are consistent with, and repro- 
duce those of [19] when T = T Ci00 . 
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FIG. 1: Scaling functions 6( pcr ) for n = 1,2, 3, oo obtained 
by extrapolating the O(e) results of [8] to d = 3, compared 
with the exact large-n result for d = 3 [5, 20]. 



Let us note that KD's two- loop results for these bound- 
ary conditions, though well-defined down to T CiOC , gave 
clear indications of existing problems. To see this, con- 
sider the scaling functions 0( por ) for n = 1,2,3, oo dis- 
played in Fig. 1, which were obtained by extrapolating 
their 0(e) results to d = 3. 

The behavior of these curves at small L differs in a 
qualitative fashion from that of the exact scaling func- 
tion for n = oo and d = 3, which follows from the ex- 
act solution of the mean spherical model under periodic 
boundary conditions [5, 21, 22]. Unlike the latter, which 
decreases monotonically to its critical value 



A (pe 



•,SM) 



-2C(3) 
5ir 



~ -0.15305 



(1.8) 



at L = 0, the former go through a minimum at small L > 
and then increase as L — ► 0. Such a minimum is neither 
expected at L > nor in conformity with the Monte Carlo 
results of Ref. [23] and announced more recent ones [24- 
26]. Note also that as n increases, the extrapolations 
actually move away from the exact n = oo curve since 
the deviations at small L get bigger. 

A second problem was pointed to by KD: Since for fi- 
nite L no phase transition takes place at T C:00 , the free 
energy per unit area must be an analytic function of 
temperature at T c oo , which imposes conditions on the 
small-L behavior of the scaling functions 0( p )(L) (which 
will be recalled in Sec. IV C). KD found their 0(e) re- 
sults to be consistent with these conditions only in the 
considered cases of non-zero-mode boundary conditions 
p = ap, (D,D), and (D, sp). In the remaining cases of 
the zero- mode boundary conditions p = per and (sp, sp), 
these conditions turned out to be violated by terms of 
first order in e. 

The results our approach yields for the scaling func- 
tions 0( per ) do better in two regards. First, the small- 
L behavior is improved inasmuch as the Casimir ampli- 
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tudes A^? cr ^ are approached in a monotonically decreas- 
ing manner as L — > 0. Second, the order of the terms 
violating the analyticity condition is increased from 0(e) 
to 0(e 3 / 2 ). In the case of sp-sp boundary conditions, 
our results raise questions whose answers might require 
a generalization of our analysis in which the surface en- 
hancement variables are allowed to vary. As we shall 
see, the one-loop expression for the scaling function of 
the inverse finite-size susceptibility becomes negative in 
a small interval of L = L/£oo when evaluated at e = 1. 
This probably simply means that this extrapolation to 
d = 3 is not sufficiently accurate. In any case, this vio- 
lation of a necessary stability condition of the disordered 
phase is a problem even for KD's original 0(e) results. 

The remainder of this paper is organized as follows. In 
the next section we specify the model utilized in our anal- 
ysis — the 4> A theory in slab geometry. We briefly recapit- 
ulate the general fluctuating Robin boundary conditions 
it involves, its renormalization, the fixed points that are 
relevant for the subsequent analysis, and the renormal- 
ization of its free energy. In Sec. Ill we first recall the 
conventional theory of the Casimir effect based on RG 
improved perturbation theory in 4 — e bulk dimensions, 
and then discuss the problems into which it runs when 
a zero mode appears in Landau theory. This is followed 
by a detailed exposition of how these problems can be 
overcome through an appropriate reformulation of field 
theory. In Sec. IV our results for the residual free ener- 
gies and their scaling functions are presented. In Sec. V 
we employ the solution of the mean spherical model un- 
der periodic boundary conditions [5, 20, 21] for d < 4 
to show that our small-e results are in conformity with 
these exact ones in the limit n — > oo. A brief summary 
and discussion of our work is given in Sec. VI. Finally, 
there are four appendixes in which technical details are 
described. 



II. CONTINUUM MODEL, BOUNDARY 
CONDITIONS, AND BACKGROUND 

A. Definition of model 

We consider a d-dimcnsional slab of finite thickness L 



occupying the region 2J = 



space. Let Xj 



[0, L] of rf-dimensional 



1, . . . , d, be Cartesian coordinates, 



with Xd = z taken along the finite direction. We write 
the position vector x = (x\, . . . , Xd) as x = (y, z), where 
y = (x\, . . . , Xd-i) is the component along the slab. 

The Hamiltonians of the (j> 4 models we are concerned 
with are sums of a bulk and a boundary term, 

H[</>]= I C m {x)dV + f C d<a (x)dA, (2.1) 

where £<s(x) and £d<n( x ) depend on <f>(x) and its deriva- 
tives. 



We either consider periodic or free boundary condi- 
tions along the z direction. In the first case, where 



4>(x + Le d ) = <j)(x) , 



(2.2) 



there is no boundary, <90J = 0, and the boundary term 
j dm ... is absent. In the case of free boundary conditions, 
the boundary (903 = 03i U 032 is the union of 03i, the 
z = plane, and 032, the z = L plane. 
The bulk density is always given by 



1 n o 



(2.3) 



where 4>{x) = (<j> a (x)) is the n-component order param- 
eter field and <fr denotes its absolute value \4>\. 

The boundary density we utilize when considering free 
boundary conditions reads 



C m [<t>] = ^4>{xf , 



(2.4) 



where c(x), the surface enhancement variable, is allowed 
to have different values on 03 1 and 032, i-e., 



c(x) 



c\ for x e 03i, 
c 2 for x e 23 2 . 



(2.5) 



B. Boundary conditions 

Using well-known arguments [12, 13], one concludes 
from the boundary terms in the classical equations of 
motion 5H = that the derivative d n 4> along the inner 
normal n on <90J satisfies 



d n (j) a (x) = Cj <f)a(x) for x € 03 j 



(2.6) 



This is a boundary condition for Landau theory, which 
holds beyond it in an operator sense (inside of averages) . 



C. Renormalization of correlation functions 

To absorb the ultraviolet (uv) singularities of the (N + 
M)-point cumulant functions 

, N M 

g[ Z] m) ,pm ( x i' • • • ' vm) = ( n ^ ( x i) n fa* ^ 

\ i=\ fe=i 

(2.7) 

involving N interior points Xj ^ <90J and M boundary 
points yk € 903 for dimensions d < 4, bulk and boundary 
counterterms are needed, which can be chosen to corre- 
spond to the reparametrizations 

4> = z l /2 <t>R , 



T - Tc.oo = 5t = Z T , 
uN d = [i e Z u u , 



(2i 
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and 



values 



SCj Cj Cgp ^Zq Cj , 

4>\d<B 



(^i) 1/2 *|* 



(2.9) 



Here 4>\q<b means <fi(yk) at a boundary point yk, and 
ji is an arbitrary momentum scale. Further, Sf is the 
deviation of f from f Cj0 o, the critical-point value of f 
of the bulk system. In a theory regularized by a large- 
momentum cutoff A, f CiOC would diverge ~ A 2 . We prefer 
to use dimensional regularization; then f c oo vanishes in 
perturbation theory. The renormalization factors Z^, Z T , 
and Z u are standard bulk quantities. The renormaliza- 
tion factors Z\ and Z c are properties of the semi-infinite 
system that results in the limit L — > oo when c\ = c 2 
with | cx | < oo. 

We choose the factor that is absorbed in the renormal- 
ized coupling constant as [27] 



N d = 



2T(3-d/2) 



(d-2){An) d / 2 



16tt 2 



1 + 



1 - C E + ln(47r) 



e + 0(e 2 ) 



(2.10) 



where Ce is the Euler-Mascheroni constant. It differs 
from the one advocated by Schloms and Dohm (see, e.g., 
[28] and its references) by a trivial factor of 2. Ours 
agrees to zeroth order in e with (2~ d 7r~ d / 2 ), the one em- 
ployed in Ref. [12] and by KD. Therefore, all of the above 
bulk and surface renormalization factors Z^,. . . ,Z\ re- 
main the same as in [8, 12] when determined by minimal 
subtraction of poles at e = 0. Explicit two-loop expres- 
sions for these functions can be found in Eqs. (3.42a-c) 
and Eqs. (3.66a,b) of [12], or in Rcfs. [17, 18]. The ad- 
vantage of our choice of Nd is to simplify the resulting 
expressions for renormalized one-loop bulk vertex func- 
tions while leaving the renormalization factors of [12] and 
KD unchanged. 

The quantity c sp is the special value of Ci correspond- 
ing to the critical enhancement of the surface interactions 
in a semi-infinite system with surface plane 25i — i.e., 
the value at which the so-called special transition occurs 
(provided the surface dimension d — 1 is large enough to 
allow long-range surface order above T C)0O ). Analogously 
to f c ,oc, it would diverge ~ A as A — > oo in a cut-off reg- 
ularized theory, but vanishes in a perturbative approach 
based on dimensional regularization [29-31]. 



D. Fixed points 

Let u* be the infrared-stable zero of the beta function 
f3 u (u) = /id^lou, where d^\o means a derivative at fixed 
bare parameters u, f, c\, and 62 of the theory. In the 
enlarged space {r, u, c\, C2} of bulk and surface variables, 
the RG yields fixed points on the hyperplane (r, u) — 
(0,u*) located at the 9 pairs (cj,^) of the fixed-point 



00 , 

0, 

—00 



(2.11) 



These values pertain to the fixed points describing the 
critical behavior at the ordinary, special, and extraordi- 
nary transitions of the semi-infinite system. Each one 
of these fixed points is specified by a pair (k\, k 2 ) with 
Ki, K2 = ord, sp, ex of the respective surface universal- 
ity classes. Universal finite-size quantities such as the 
Casimir amplitudes and the scaling functions Q^ p \ 
~(p) generally are different, depending on the basin of 
attraction of the fixed point (ki,K2) to which they be- 
long. Recall that for kj = ord the cumulants (2.7) satisfy 
the Dirichlet boundary condition lim X( .^<8 j G^ ) = 0. 
Thus the universality classes (ord, ord), (ord, k), and 
(k, ord) with n ^ ord can equivalently be labeled as 
(D,D), (D,k), and (k,D), respectively. We shall con- 
tinue to employ this convention. 



E. Renormalization of free energy 

The counterterms implied by the reparametriza- 
tions (2.8) and (2.9) are sufficient to absorb the uv singu- 
larities of the cumulants (2.7). However, the free energy 
requires additional additive counterterms [8, 12]. They 
can be chosen to be independent of L [12, 14]. We there- 
fore add to the Hamiltonian defined by Eqs. (2.1)-(2.5) 
a contribution 



*add 



= f Cx{T,u)dV + y] 



Cd<a(f,u,Cj)dA , 



where is a polynomial in f of degree 2, 



(2.12) 
Cggj is a 



polynomial of degree one in f and degree 3 in Cj, whose 
coefficients depend on u, but neither on L nor on the 
position x. The coefficients (power series in u) are fixed 
as follows [32]. Let T^p/(f) denote the Taylor series 
expansions of the function / to second order in f (4th 
order in f 1 / 2 ), and 



Tgg(r,c) = 



E 

0<2j+fc<3 



1 di+ k g 



>\k\ d¥dc k 



NP 



x (f-r N p)^(c-c N p) fe . (2.13) 

be the corresponding expansion of g in f and c to orders 
j and k with 2 j + k < 3 about the normalization point 
(f, c) = (tnp,cnp) [33]. We choose 

,2 Z T , 



T"NP = T\ T= i = T c . 



o +M 



c NP =c\ c=1 = c sp + {iZ c , (2.14) 

and define the dimensionless renormalized bulk free en- 
ergy density f btR by 

h,R{T,u) = fb{r,u) 



(2.15) 
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The excess surface free energy density f s (f, u, ci, c 2 ) 
of the infinitely thick film is a sum of contributions 
f s (f,u,ci) associated with the respective semi-infinite 
systems bounded on one side by <Bj-; i.e., 

/ s (f,u,ci,c 2 ) = / s (f,u,ci) + f s (f,u,c 2 ) . (2.16) 

We define the dimensionless renormalized analogs of the 
latter by 

f s , R (T,u,c) = / s (f,w,c) - T^pf s (f,U,C) . 

(2-17) 

By construction, these renormalized bulk and surface 
free energy densities satisfy the normalization conditions 



fb, R \ 



NP 



/ 6ifl (l,u) = = 



df b ,i 



dr 



d 2 f, 



b,R 



NP 



Or 2 



NP 

(2.18) 



and 

d j+k f s ,R 



NP 



(1,«,1) = 0, 0<2j + fc<3, 



3 T j ,9 c fc 

(2.19) 

respectively. The renormalization functions C<jj and Cgaj 
are fixed by these requirements. 

The renormalized excess surface free energy density 
/ s ,ij(r, u, Ci, c 2 ) one obtains from the action 7i + ^4 a dd 
upon insertion of the reparametrizations is uv finite. 
Since Caj and C^qj are independent of L, the subtrac- 
tions they provide cancel in the residual free energy 
f les (L;f,u,Ci,c 2 ) of the film. Accordingly, its dimen- 
sionless renormalized counterpart 

fr CS .R(pL;T,u,c 1 ,c 2 ) = ^ _(d_1) / ros (L;f,w, ci,c 2 ) 

(2.20) 

satisfies a homogeneous RG equation, whereas both f b<R 
and / Sj ij satisfy inhomogeneous ones. 

Following the notation conventions of Ref. [12], we in- 
troduce the beta function f} u = ^id^ou, the RG functions 
Vk = Atd^lornZjt, k = <f), t, m,c, 1, and the operator 

2 

^ = M^+/?A-(2+r ?r )9 r -(l+r ?c )^c j a c . . (2.21) 

Then the RG equation of / ros ,ii can be written as 

[X> M + {d- l)]f, cs . R (pL; t, u, cud) = . (2.22) 

Note that the RG functions (3 U and r\ K are either bulk 
quantities (such as /3 u (e, u), rj T (u)) or properties of semi- 
infinite systems such as ij c (u). In accordance with 
Ref. [12], we have chosen them independent of Cj and 
t (fixing them by minimal subtraction of poles in e). 
Explicit two-loop expressions for these functions can be 
found in Eqs. (3.75a)-(3.76b) of this reference. 

The RG equation (2.22) can be solved in a standard 
fashion by means of characteristics. Upon setting [i = 1 
and choosing the scale parameter £ of the transforma- 
tion \i — > /i£ equal to l/£oo, the inverse bulk correlation 
length, we see that the residual free energy density, on 



sufficiently long length scales, takes the finite-size scaling 
form 

f les , R (L\ r, u, d,d) « Q(L/U,d C,c 2 C) ■ 

(2.23) 

Here $ is the surface crossover exponent of the special 
transition. The scaling function O is universal up to the 
nonuniversal amplitude of ^ and the nonuniversal met- 
ric factor associated with Ci and c 2 [34] ; it is given by 

9(L; ci, c 2 ) = L^ 1 / res , fl (L; l,u*, ci, c 2 ) . (2.24) 

From it the functions 0( p )(L) with p = (sp,sp), 
(D,D), and (sp, D) follow by setting (ci,c 2 ) to the re- 
spective fixed-point values (0,0), (00,00), and (0, 00). 
For example, 

e (sp,sp) (L) = \_d-i f Tes R (\_- i j u *,0, 0) . (2.25) 

For reasons explained in the introduction, we shall 
mainly be concerned with the cases of periodic and 
(sp,sp) boundary conditions. 



III. REVISED FIELD THEORY APPROACH 
A. Infrared problems due to zero modes 

We now turn to the problem of computing the scaling 
functions 0^ and for p = per, (sp, sp) by means 
of RG-improved perturbation theory. Only the case T > 
T c ,oo will be considered. 

The free propagator can be written as 

^ (z\m)(m\z') 



G^(x;x'\f) = f " ]T 



P 2 + k 2 m + t 



,ip-(y-y') 



where 



(3.1) 



is a convenient shorthand for a normalized d — 1 dimen- 
sional momentum integral. Further, |to) = \m)^ are 
eigenstates given by 



{z\m} 



(per) _ 



cxp(ifc m z) 



and 



(z\m) 



(sp,sp) 



J- J' 



k m — ^^- W , m e Z , (3.3) 



for to = 



VL \y/2 cos(fc m z) for to e N , 
k m = mn/L, m = 0, 1, . . . , 00 , (3.4) 



respectively. For cither boundary condition the mode 
with to = and p = becomes massless at r = 0. 

In their calculation of A^f* directly at T c oo , KD there- 
fore subtracted the contribution from the to = mode 
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to avoid infrared problems, the rationale being that the 
subtracted one-loop contribution is formally independent 
of L so that it does not contribute to the Casimir force. 
Computing the one- and two- loop graphs at r > [8, 19], 
one finds that the contributions from the ko = modes 
vary as a positive power of r and hence vanish as r — * 0. 
However, at the three-loop level this is no longer the case 
because one encounters infrared divergent contributions 
of the form depicted in Fig. 2. Thus conventional RG- 



FIG. 2: (color online) Infrared divergent contribution to the 
free energy. The dashed full blue lines represent the k m 7^ 
part of the free propagator (3.1); the dotted red lines denote 
its k m — part. The blue subgraphs approach a finite In- 
dependent limit as r — > 0; the red dashed subgraphs varies as 
a negative power of r and hence is infrared singular [19]. 

improved perturbation theory is ill-defined at T c oo . 

The origin of this problem is that Landau theory yields 
sharp transitions for both the bulk and the film system at 
the same critical value f = 0. It is thus of a similar kind 
as encountered in the study of finite-size effects of systems 
that are finite in all, or in all but one, direction under 
periodic boundary conditions [35, 36]. As discussed in 
Ref. [19], the remedy is to separate the fco = mode and 
construct an effective field theory for the fco = part of 
the order parameter. 



B. Construction of effective zero- mode action 

To this end we write 

<j>{ x ) = Y, <M2>) (*M = L ~ 1/2 rtv) +1 Kv> z ) ' ( 3 - 5 ) 

m 

decomposing the order parameter into its component 
along 4>o(y) = <p{y) and a remaining k m ^ contri- 
bution "0(y, z) with 



/ dzip(y,z) = Q 
Jo 



(3.6) 



Tracing out ip defines us a (d — l)-dimensional effective 
field theory with the Hamiltonian 



1~(e«[tp} = -InTr 



-HlL- x / 2 <p+4,] 



H[L-^ 2 ip] - ln{e 



-Hint[<P,i>]\ 



k B T 

Here , defined by 

exp(-iVfc B T) = Tr^cxp (-H[ip]) , (3 
is the free energy due to the k m ^ modes. Further, 



(3.7) 



H ia t[v>,il>] = / dV 



„2 „/,2 



( V -ip)i> 2 (3.9) 



is the interaction part, and 
d-i 

H[L- X l 2 v\ = I dA" 



V (&p_ y 
2 ^ KdxJ 

3=1 3 



2^ 



AIL 
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(3-10) 

Computing the last term on the right-hand side of 
Eq. (3.7) in a loop expansion gives 



with 



F, 



■4, 



knT 



Tr In 



■H[L-V 2 <p}+nW[<p] 



(3.11) 



1 + ^G^(6 af3 ip 2 +2<p a tp l3 ) 



I I LP . LD 

(3.12) 

where the dashed blue lines (color online) represent free 
■0-propagators 

,(d-i) 

■ e" 



(Z\m)(m\z') ip .(y-y>) 



P 



U2 



and the red bars indicate Lp legs. 
Writing 

5 2 UM 



(3.13) 



<W g^{y - y') 



5(p a (y) 5(pp(y') 



<p=0 



= 6 aP [(-d 2 + r)s(y-y')-^(y-y')] .(3.14) 

where d 2 means the Laplacian m '^■ d ~ 1 , and 



(v-i 

!oti,...,a k Vtflj 



,Vk) 



8 k n cS {ip\ 



8<p ai (yi) ■ ■■8 ( p ak (yk) 



ip=0 

(3.15) 

we introduce the propagator g v associated with the if 2 
term of 7i c ff [v 3 ] an< l t ne corresponding self-energy a v as 
well as the vertices 7'*) of the effective action H e s[<p]. 
Though not indicated here, all these quantities depend 
on L and the boundary condition p. 

The first graph in the second line of Eq. (3.12) is the 
one-loop contribution to a v . It is local in y-space. As 
can be seen from Fig. 3, both local and nonlocal contribu- 
tions appear beyond one-loop order. The second graph 



ip. 



FIG. 3: (color online) Two-loop contributions to oy. The left 
graph is local, the one on the right-hand side is nonlocal. 

in the lower line of Eq. (3.12) is the nonlocal one-loop 
contribution to j( 4 \ Evidently, vertices 7'*' of arbitrary 
even order k are generated through the coupling to the 
k m 7^ modes. 
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C. RG-improved perturbation theory 



where 



Now suppose the bulk critical point is approached so 
that £oo becomes large. Then the vertices 7^ cannot be 
computed by perturbation theory below the upper criti- 
cal dimension d* = 4. However, for arbitrary small r > 
we can employ the RG to map to a system with a min- 
imal length scale on the order of £oo, and then employ 
perturbation theory. The vertex functions 7^ are ex- 
pected to decay as a function of the relative differences 
Vij = \y%-yj\ on the scale of . 

The renormalized counterparts 7^ of these vertices 
satisfy the RG equations [37] 



N 



(fe) 



. 



(3.16) 



Solving them in a standard fashion, one finds that the 
Fourier transforms 7^ (27r) d_1 5Q2j =1 Pj) of these func- 
tions on sufficiently large length scales take the scaling 
forms 

iT\{Pi})~^^-w-^ k x^{{ Pi ^y,L/^), 

(3.17) 

where rj is a standard bulk critical exponent, while ^ is 
the second-moment bulk correlation length. The latter is 
defined in the conventional manner in terms of the bulk 

~ (2) ~ (2) A 

vertex function 1^ (q) = 1/G h (q) of the <j) theory in d- 
dimcnsional momentum space or its position-space back 
transform (x) via 



d 



£2 In if ^ 



1 fd*xx*Gg\x) 



q= o 2d jd d xG^\x) 

(3.18) 

Let us verify explicitly to first order in u* = 0(e) that 
RG-improved perturbation theory yields such scaling be- 
havior. Consider, for example, 7^ = The first 
graph in the second line of Eq. (3.12) is the 0{u) contri- 
bution to <r v (p). Introducing 



(L;f) = / 
Jo 

we have 



dz 



G[%(x;x\f) , . 7 = 1,2, (3.19) 



-u^I^(L;f) + 0(u 2 ). 



(3.20) 



2N d 
4-d 



.(4 7r )-d/2 r ( 1 -d/2). (3.23) 



Here Qd,2(r) is a special one of the functions defined by 



Qd,*(r) 



E 



(d-1) 



(d) 
q=(p,k) 



^a-2 



q 1 + r 



(3.24) 



where p and k are the (d— l)-parallel and one-dimensional 
perpendicular components of the wave vector q — (p, k). 
The properties of these functions are analyzed and dis- 
cussed in Appendix D, where we compute them for the 
required parameter values of d and a. Plots of the func- 
tions Qdfiir) with d — 4 and 6 are displayed in Fig. 9 
(Appendix D). 

To facilitate subsequent comparisons with KD's re- 
sults, let us note how the Qd,2{r) are related to the func- 
tions 



9a,b{z) 



a J 1 



dt 



(t 2 - If \n\t 2 - 1) 



a 2zt 



1 



(3.25) 



utilized by these authors. As shown in Appendix C, one 
has 

2 1 - d 7T ( 1 - d )/ 2 r d / 2 
Q d .2(r) = r[(d _ 3)/2)] ^,o(VF/2) . (3.26) 

Using the results (3.21) and (3.22), and expressing 
9<f ii = 9<p 1 m terms of the renormalized variables 
r and u, one finds that the pole <~ e _1 cancels. The re- 
sulting renormalized expression is easily evaluated at the 
fixed-point value u = u* . It conforms with the scaling 
form 

[g^ R {p,r,L)}-^ fjP L^ 2 X^\p^L/U) . (3.27) 
and yields for the scaling functions the e expansions 



Af c %,L) = (p 2 + l)L 



[2ttL 



and 



X. 



(sp,sp) 



(P,L) 



n + 8 

+ 16^ 2 g 4 , 2 (L 2 )/L 2 ] +Q(e 2 ) (3.28) 



( P " + l)L" + 2±|«[.L 

+ ^ 2 Q 4 , 2 (4L 2 )/L 2 ] +0(e 2 ) . (3.29) 



The integrals Ij (L;f) are computed in Appendix A. 
The results for ij pcr ) and ij sp ' sp ) are 

j(pcr) (£; f ) = Ad ~ X f (d ~ 3)/2 - A d T (d - 2)/2 + 2 ® d ^ L2 ) 

(3.21) 

and 



/(^.»p)(L;f)= j(P er )(2L;f) , 



(3.22) 



A few comments are in order here. 

(i) The above results imply that [fl , £f/j(0, r, docs 
not vanish at r = when L < 00. Using the fact that 
rj = 0(e 2 ) and the small-r behavior of implied by 

Eq. (A15) yields 

[ff^ sp) (0,0,L)*]- 1 = [4^(0 ) 0,L)*]- 1 +O(6 2 ) 



_2 C(2) 
8 L 2 



+ C>(e 2 ), (3.30) 
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where the asterisk indicates evaluation at u = u* . 

The physical meaning of this result is obvious. The 
coupling of the k = mode to the k m ^ modes has 
produced an L-dcpendent shift of the temperature at 
which (p becomes critical, making ip noncritical at T C:00 
when L < oo. 

(ii) Verifying the scaling form (3.27) to higher orders 
in e and the appearance of a nontrivial exponent n by 
extending RG-improved perturbation theory to 0[(tt*) 2 ] 
or higher is in principle straightforward. 

(iii) It is instructive to see what our procedure yields 
for boundary conditions such as p = (D, D), (D, sp), and 
(ap) where Landau theory does not involve a zero-mode 
at T CiOC . In those cases we have if = and <p = ip. Ac- 
cordingly, Eq. (3.11) simply yields F/fc B T = F^/ksT for 
the reduced free energy. It is therefore clear that for those 
non-zero-mode boundary conditions conventional expan- 
sions in integer powers of e will result for the Casimir 
force, the scaling functions and Q( p \ and similar 
quantities for T > T C:00 , which must be in accordance 
with KD's results to 0(e). 

(iv) The procedure utilized above of constructing the 
action of an effective lower-dimensional field theory by in- 
tegrating out modes via RG-improved perturbation the- 
ory that do not become critical for T = T CjOC at zero- 
loop order is similar to the one employed in the study 
of static and dynamic finite-size effects in systems that 
are finite in all, or in all but one, directions [35, 38- 
41]. In the latter cases one arrives for small deviations 
e = d* — d > from the upper critical dimension d* = 4 
at expansions in powers of e 1 / 2 and e 1 / 3 , respectively. 
The main difference between these cases and ours is that 
a sharp transition to a low-temperature phase with long- 
range order is ruled out for the former because they in- 
volve systems of finite extent along d or d — 1 Carte- 
sian axes (and the presumed short-range interactions). 
By contrast, in the case of the slab geometry consid- 
ered here, such a sharp transition should occur for finite 
thickness L at a shifted temperature T c ,l < T c ,oo when- 
ever d—1, the effective dimensionality, is sufficiently large 
for such a long-range ordered low-temperature phase to 
occur. (Evidently d—1 must exceed d*(n), the lower 
critical dimension, which is d*(l) = 1 in the Ising case 
n = 1, and d*(n > 1) = 2, depending on whether a 
discrete Z 2 or continuous 0(n) symmetry gets sponta- 
neously broken.) When no sharp transition is possible, 
one expects a rounded one at a shifted pseudo-critical 
temperature (see, e.g., Refs. [42, 43]). The case of d = 3 
and n = 2, corresponding to an AY~-model on a slab 
or liquid 4 He film below the bulk A-line T\, is excep- 
tional in that a transition of Kostcrlitz-Thouless type to 
a low-temperature phase with quasi-long-range order is 
expected to occur for finite L. 

(v) That the coupling of the fco = mode ip to the 
k 7^ modes ip produces an L-dependent mass gap for 
g^ 1 is crucial for making RG-improved perturbation the- 
ory well-defined at T Cj00 . However, it must be empha- 
sized that such a perturbative approach using u* = 0(e) 



as expansion parameter by itself must not be expected 
to give a proper description of the (d—1 dimensional) 
critical behavior at T c j}- One way to see this is to note 
that the bare ip A coupling constant appearing in H c s [<fi] 
is ii/L. To make it dimensionless we must multiply by 
the (5 — d)th power of a length. An appropriate one is 
£l, the finite-size analog of £oo, defined by 



f 2 



p=0 

2 



1 /^-WMy)-v(0)) c 



2(d-l) / d^y (<p(y) ■ <p(0)Y 



(3.31) 



where f w (p) denotes the full ipip vertex function in the 
space of (d — l)-dimensional momenta p. 

The appropriate dimensionless coupling constant 
therefore is ^~ d u/L, which diverges as £l — ► oo when- 
ever d < 5. In accordance with general expectations we 
thus see that the appropriate smallness parameter for an- 
alyzing the d—1 dimensional critical behavior at T c .l by 
means of a dimensionality expansion is b — d rather than 
e. Constructing a RG approach that is reliable both at 
T c oo and T c ^ and capable of describing the crossover 
from d to d — 1 dimensional critical behavior is a non- 
trivial problem, which has so far not been solved in a 
satisfactory fashion and is beyond the scope of this pa- 
per. 



IV. CALCULATION OF FREE ENERGIES AND 
SCALING FUNCTIONS 

According to Eq. (3.11), the reduced bare free energy 
density per unit area 



n 



F 



lim 
a^oo Ak B T 



is a sum 



/r(r) = 4 P) (£;r) + /^(L;f) 



(p)i 



(4.1) 



(4.2) 



of a contribution f^'(L; f) from the k m ^ modes and 
a remainder, which we denote as f!f\ We first consider 
the non-zero mode contribution . 

A. Non-zero mode contribution to the free energy 

A standard loop expansion yields 
f { /\L;f) = 4^ ] (L;f) + /^ [ ) 2] (L;f) + 0(3-loops) (4.3) 
with 

r(d-l) 



= 4^ ] (L;0) + L->)(L;f) (4.4) 



9 



and 



where 



uL 



n(n + 2) (p) 



4! 



I^(L;f) 



(4.5) 



J (p) (L;f)= / l[ p) {L;t)dt. (4.6) 
Jo 

The f = contributions /^, P [\] (£; 0) arc computed in Ap- 
pendix B. The results arc in accordance with those of 
KD. Expressed in terms of the familiar one-loop values 

Ag~} = 2 d Agjif > = -nn-^ T(d/2) ((d) (4.7) 

of the Casimir amplitudes, they can be written as 

/S]^°) = 4^+ i " (<< " 1)A c1i]- (4-8) 

Here are the cut-off and L dependent quantities de- 
fined by Eqs. (B3)-(B5); they vanish in dimensional reg- 
ularization. 

The integrals I { 2 p) and J(p) are worked out in Ap- 
pendix A. The results are given in Eqs. (A4), (A9), 
(A16), and (A17). Inserting them into Eqs. (4.4) and 
(4.5) gives 



A d . 



I f(d-l)/2 _ L ^± j.d/2 
1 d 



47rQ d+2j2 (fL 2 ) 



f^(L;f) = uL 



-L d+1 



n(n + 2) 
4! 



(4.9) 



A 



L 



and 



f (sp,S P ) ,j „x _ f (sp,S P ) 



1 Arf 



2d- 1 



2Q d ,2(rL 2 ) 



-±f(d-l)/2 



(4.10) 



d T 2 d iL d + 1 



, (4.11) 



f (SP,S P ) / j . o x 



uL- 



i(n + 2) 



A 



d-l f (d-3)/2 



2L 



■j .(d-2)/2 , Q^2(4t£ 2 
AdT + 2 d ^fL d 



+ 



+ 



B d 



A 



+ 



2L 



2d-3 T 



<5 2 d-3,2(4fL 2 



22d-4^£2d-3 



(4.12) 



respectively, where /jfg and B d are constants defined in 
Eqs. (B3) and (A10). ' 



The contribut 



ibutions —nLA d T d / 2 /d to /^fm have sim- 
ple poles at d — 4, which get cancelled upon renormal- 
ization by the additive bulk counterterm cx r 2 implied 
by the subtraction (2.15). The two- loop terms f^ 2 ] m " 
volve uv singular bulk terms linear in A d whose poles 
at e = get cancelled by the 0(u) contribution to the 
counterterm (Z^ Z T — 1) /i 2 t <j) 2 R /2. That no pole-term 
singularities located at the boundary planes Q3i and 2$ 2 
appear at e = in /^j^ and is because both 

renormalized enhancement variables c\ and c 2 are zero. 

Since our main interest is in the renormalized residual 
free energy f^ R , we can avoid dealing with additive 
countcrterms by focusing directly on its calculation. To 



(p) 

*/>,res,.R> 



we 



determine its non-zero mode contributions 
must subtract from the sums of the above one- and two- 
loop terms the bulk and surface contributions and express 
the difference in terms of the renormalized variables r 
and u: 



= 4 p) (L;f,&)-L/ 6 (p) (f,fi)-45(f,fi). (4.13) 

From the results (4.9)-(4.12) one easily reads off the 
(p-independent) bulk terms 



fb — fb,o 



nA. 



.r d / 2 +u n{n + 2) A 2 r d - 2 + 0(u 2 ) 



d 4! 
as well as the p-dependent surface terms 

p(per) _ f(pcr) , n A d _i ,°(d-l)/2 



(4.14) 



_ „ n(n + 2) fd _ 5/2 + Q(fl2) (4 l5) 



and 



f (sp,sp) _ 1 ,(pcr) „ n(n + 2) °d-5/2 j_n(n°, 2 \ 

Jip.s ~2^' s 4! B d A 2d-3T +U{u ). 

(4.16) 

No confusion should arise from the fact that /^^^ does 

not vanish. As is easily checked, and our results for flf 01 " 1 
to be given below will show, this term cancels exactly 
with the surface contribution to f!f CI \ as it must. Of 
course, such cancellations are neither expected nor occur 
for p = (sp,sp) and other boundary conditions. 

With the aid of the property 



d_ 

dr 



Qd+2,2( r ) 



QdMr) 

471T 



(4.17) 



derived in appendix B of Rcf. [21], the calculation of 
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fib res r becomes straightforward, giving 



„(pcr) 
J ib,res,R 

n 



tL 2 ) u n + 2 



■(jjlL) 



d+1 



\xL 4! 



N d 



Ad i r (d-3)/2 + ^Qd,2{p 2 ' rL2 ) 



d-1 



r-^~l8Q d . 2 (fi 2 rL 2 ) 



d-i 



■QmL) 
\+0(u 2 ) 



(4.18) 



and 

As 

' ib,res,R 



f(sp,sp) 
Jib 



2tt Q d+2 a{^ 2 TL 2 ) u n + 2 



fi d 1 n 



T(2 f iL) d + 1 



+ 



fiL 4! 



\N d 

Bd 
N d 



Ad-i (d -3)/2 Qd,a(4M 2 TL 2 ) 
2 T r(2 M L) d - 1 



^2d-4 d _3 , Q2d-3,2(4M 2 Ti 2 
T 4 



2 ' T(2/zL) 2d - 4 
r- £ / 2 -12g d , 2 (4 M 2 ri 2 



d-2 




(2/iL) 



B. Remaining free energy terms 

We next turn to the computation of f^. For u = 
the Hamiltonian H c s[(p] describes a free field theory 
whose two-point function is the familiar Gaussian bulk 
propagator 



(d-i) 



{p 2 +r)- 1 e lp - y 



2 Wd-3)/4 



jy^K^y^yVf) (4.20) 



(27r)( d -!)/ 

in d — 1 dimensions. As we have seen, using this as free 
propagator in a Feynman graph expansion would lead 
to Feynman integrals that are infrared divergent at T c oo 
and make the expansion ill-defined beyond two-loop or- 
der. This suggests to work with a free propagator whose 
mass parameter, firstly, remains positive for T > T c oo 
when L < oo, and secondly, has a well-defined physical 
meaning beyond perturbation theory. A natural candi- 
date that has these properties is the inverse finite-size 
susceptibility = rt, defined by 



(r L ) 1 <W = XL5 af3 = / d d 1 y{i Pa (y)Lp p {Q)) 



We therefore use 



GJy) = G<£ 



(4.21) 
(4.22) 



as free propagator. A tacit assumption underlying our 
calculation is that the disordered phase is the correct 
reference state to expand about for the parameter values 
of L and t > considered. Since the transition tem- 
perature T C] l for finite L is expected to be lower than 
the bulk critical temperature T Cj00 (in those cases of d 
and n for which a sharp transition occurs when L < oo), 
this is physically reasonable. However, there is no a priori 
guarantee that extrapolations to d = 3 of results based on 
RG-improved perturbation theory will fulfill all necessary 
requirements. In particular, we should check whether the 
so-obtained approximate inverse finite-size susceptibili- 
ties remain positive when L < oo. This issue may 
be expected to be more delicate for p = (sp, sp) than for 
periodic boundary conditions. The reason is that sp-sp 
boundary conditions are associated with a multicritical 
point of the surface phase diagram (located at c = t = 0) 
at which the line of surface transition T CjS (c) meets the 
bulk critical line (whose sections with c > and c < 
form the lines of ordinary and extraordinary transitions, 
respectively) [12, 44]. For finite L, one expects shifts of 
this multicritical point and the phase boundaries. To ac- 
count for these shifts one would have to vary the surface 
enhancement variables Cj as well, giving up the restric- 
tion c\ = C2 = 0. This is a difficult problem and beyond 
the scope of the present investigation. 

Let us represent the propagator (4.22) by a red line, 
the effective two-point vertex f — fx — a v by a red dot 
with two legs, and the effective fc-point vertices with 
k > 2 by red dots with k legs. Then the Feynman graph 

expansion of f!f^ (L) becomes 



-fi P \L) 




(4.23) 



The first graph on the right-hand side is given by 

- o = n f *>-° + \ r g -" i} ( ° ^ ^ (4 - 24) 



with 



fvfi - 



(d-i) 



hip 2 



(4.25) 



Our results (3.14) and (3.20) for g 1 and a v imply that 



r-r L -a v = H±* | G^ 1 ) (0\f L ) + 0(u 2 ) 



n + 2 A d _iu (d_ 3 )/2 
6 L T L 



r[ d - 3)/2 + 0(u 2 ) . (4.26) 



Using this in conjunction with the fact that the effective 
four-point vertex, to first order in u, is a local (p 4 cou- 
pling with interaction constant u/L, one finds that the 
contributions from the other two graphs can be written 
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as 



where 



u n(n + 2 
L 



4! 



G^(0\r L )\ +0(u 2 ) 
lQ+0(£ 2 ). (4.27) 



Upon inserting the y — > limit of the free Gaussian 
propagator (4.20) into Eqs. (4.24) and (4.27), the re- 
quired integrals can be performed to obtain 



A, 



d-l o(d-l)/2 



u n(n + 2) a2 



4! 



£.d-Z 



d-l 1 L 



(4.28) 



Just as the constants introduced in Appendix B, 
f v fi involves uv divergent contributions which arc elimi- 
nated in the renormalized theory by the additive renor- 
malization of the free energy. Furthermore, it should be 



remembered that ft ■ 
condition p. We have 



depends on the boundary 



s(p cr ) 



ulAA {d - 2)/2 - 



ig{L 2 f) 
fL d 



and 



+ 0{u 2 ) 
(4.29) 



°( S P. S P) 



r(2L) d 



A 



2L 



(4.30) 



respectively. 

The 0(u) contributions oc Ad — Nd/e of have uv 
poles at e = 0. These are cured by the bulk counterterm 
(Zq Z T — 1) /x 2 t Jj, 4>%/2. For the resulting renormalized 

dimensionless inverse susceptibilities v ® = Z^ rff* / H 2 
one obtains 



„(p cr ) _ 



G 



2 u 



e/2 _ 1 



2Q d , 2 (v 2 TL 2 ) 



(li 2 L 2 r) d / 2 N d 



6/2 

+ 0{u 2 ) 



(4.31) 



and 



„0p> s p) 



2 u 

~ W2 



r e/2 



Ai-1 



2Q d , 2 (4 f i 2 rL 2 



{4p 2 L 2 r) d / 2 N d 



e/2 2pLr 1 / 2 N d 
+ 0{u 2 ). (4.32) 



In Sec. IV D we will verify that these results comply 
with the scaling form 



r[ p) = r oa RW(L/&o) 



(4.33) 



(4.34) 



is the inverse bulk susceptibility, and try to employ them 
to determine the scaling functions by means of the 
e expansion. 

Returning to the calculation of free energies, we now 

subtract from f!p in Eq. (4.28) the surface contribution 

f(p)i 



ftp \L=oo,TL-=raa ^° obtain the associated contributions 
/yfres to the residual free energies. Expressing the result 
in terms of renormalized quantities then yields 



f(p) 



A 



d-l 



d-l 



r (d-l)/2 _ (d-l)/2 



u n + 2 A 2 



r d L ~ 3 + 0{u 2 ) (4.35) 



4! N d 
for their renormalized analogs. 

C. General properties of the scaling functions 

Before we embark on the calculation of the scaling 
functions R( p )(L) and 0( p )(L) of the inverse finite-size 
susceptibility and the residual free energy, it will be help- 
ful to discuss some general properties they should have. 

In the limits L — > oo and L — > 0, Eq. (4.33) must yield 
the correct bulk behavior and finite positive finite-size 
susceptibility, respectively. This implies 



r(p)(L) 



1 



for L 



p$ l r '- 2 with p$ > , for L -> . 

(4.36) 

Turning to the free-energy scaling functions 0( p )(L), 
let us first consider their limiting behavior as L — > 0. This 
must comply with the requirement that the finite-size free 
energy be analytic in T at the bulk critical temperature 
when L < oo. As explained by KD, this translates into 
the limiting form 



0(p)(L) 



a b+ L d 



,(p) I d-l 



l- *o a(l-a)(2-a) a s (l - a s )(2 - a s ) 

(4.37) 



k=l 



where a s — a + v is a familiar surface critical index (of 
the surface excess specific heat [12, 44]). Further, at + is 
a universal number whose e expansion 



32^ 



1 



ln(47r) - C E + 3 • 



0{e 2 ) 



(4.38) 

may be gleaned from equation (8.12) of Ref. [8]. The plus 
a s+ , and as usual indicate that these 



signs at a&+ 

numbers pertain to the limit r 



0- 
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The first two terms on the right-hand side of Eq. (4.37) 
remove the singularities of the subtracted bulk and sur- 
face contributions to f^ R ; the remaining power series 
involves integer powers of t oc (T — T Cj00 )/T Ct00 . Note 
that neither nonlinear contributions to the temperature 
scaling field have been taken into account nor those of 
irrelevant bulk and surface scaling fields. Both sources 
would entail corrections to the leading thermal singular- 
ities of the bulk and surface free energies. The implied 
additional terms nonanalytic in temperature would have 
to be removed as well in the finite-size free energy and 
hence entail further nonanalytic contributions to the lim- 
iting small-L form (4.37). 

The absence of boundaries in the case of periodic 
boundary conditions implies that the surface amplitudes 

af^ are exactly zero. For the other case of interest, 
p = (sp, sp), one has 



(sp.sp) _ 



1287T 



1 + e 



2 + hi7r - C E + ■ 



n + ; 



+ G(e 2 ) 



(4.39) 

according to KD's equations (E6) and (E9). 

We next turn to a discussion of the limiting forms of 
the functions <d^(L) as L — > oo. Since we have chosen 
periodic boundary conditions along all d — 1 parallel di- 
rections yj, no edge contributions ~ L d ~ 2 to the total 
free energy are expected. Accordingly the residual free 
energy should decay exponentially as L = L/£oc - * 00 • 
The asymptotic behavior should simply follow from per- 
turbation theory. 

To become more precise, it is useful to recall the 
representations (see, e.g., equations (4.2) and (4.12) of 
Ref. [12]) 

OO 

G { r\x 12 \f)= G^(x 12 +jLe z \r) (4.40) 

j=—oo 

and 

oo 

G^ sp \ Xl ,x 2 \r) = ]T [G^(x 12 + 2jLe z \r) 

j=—oo 

+ G^(x 12 + 2(jL + z 2 )e z \f)] 

(4.41) 

of the free propagators in terms of the bulk propagator 
g£\ where x 12 = x\ - x 2 = (yi,zi) - {y 2 , z 2 ). 

The j = terms G^ 1 (xi 2 \f) yield the bulk contribu- 
tions of G L p) . The j = term G { £ {x l2 + 2z 2 \f) and 
the j = -1 term G ( £(x 12 + 2{z 2 - L)\t) in Eq. (4.41) 

represent surface contributions. Since G < £ > (x\f) decays 
exponentially as \x\ — > oo it is clear that of the remaining 
terms those involving spatial differences that are con- 
strained by the smallest lower bounds will govern the 
limiting large-L behavior of the functions Q^ p \ In the 
case of periodic boundary conditions, this applies to the 
j = ±1 terms, which involve position vectors of lengths 



> L. Hence 9( per )(L) must vary as <~ e L in the limit 
L — > oo, up to powers of L. 

On the other hand, for p = (sp, sp), there are four 
contributions involving position vectors constrained by 
the lower-distance bound 2L which govern the large-L 
behavior. Thus 6( sp ' sp )(L) must decay ~ e~ 2L , up to 
powers of L. 

To elaborate on these arguments, one can employ the 
above expressions (4.40) and (4.41) for the free propaga- 
tors in perturbation theory, dropping all of their sum- 
mands that do not contribute to the leading large-L 
behavior. In the case of the one-loop integrals it is 
again convenient to first determine the large-L forms of 
their f-derivatives and then integrate with respect to 
f . However, from our perturbative results gathered in 
Eqs. (4.9)-(4.12), (4.24), (4.27), and (4.35), the one- and 
two-loop Feynman integrals with all contributions to the 
free propagators included can be inferred. Thus no re- 
newed calculation is necessary. To determine the large-L 
behavior of the 0( p ) we must merely replace the func- 
tions Qd, 2 and Qd+ 2 , 2 by their asymptotic forms (D4) 
given in Appendix D. This yields 



e (per) (L) 



l^oo (2 7 r)( d - 1 )/ 2 



L (d-l)/2 e -L [1 + 



and 



0(u*)] 
(4.42) 



@(sp,sp)^ 



2 d JT (d-l)/ 2 



|_(d-l)/2 e -2L 



[l + 0(«*)]. 

(4.43) 

Finally, let us briefly recall what can be said about 
the behavior of the t < analogs of the scaling func- 
tions 9( p \ which we denote as 6^(L), at the transi- 
tion temperature T c l of the film in those cases where a 
sharp transition to a long-range ordered phase is possi- 
ble for finite L, such as in the Ising case n = 1 for bulk 
dimension d — 3. As a function of the temperature de- 
viation tL = (T — T c x)/T c ,l, the excess free energy per 
cross-sectional area A must have a contribution that be- 
haves as <~ t ] ^ ad ~ 1 as tL — * 0, where ctd-i is the specific 
heat exponent for bulk dimension d — 1. The transition 
point translates into a nonzero value Lo at which the func- 
tions eL p) ( L ) behave in a nonanalytic fashion. Standard 
matching of the temperature singularities then yields the 
behavior 



(4.44) 



where v = v^, as before, is the correlation-length expo- 
nent of the d-dimensional bulk system. 

A well-known consequence is that the critical- 
temperature shift varies as [42, 43] 



T C , L )/T C . 



(4.45) 



This conclusion that the shift exponent is given by \ jv 
is more or less automatic when the finite-size scaling 
form (1.5) of the residual free energy applies and hence 
is in complete accordance with our theory. 
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D. Scaling functions of inverse finite-size 
susceptibilities 

We proceed by combining our perturbative results of 
Sec. IV A and IV B with the RG to compute the desired 
scaling functions, beginning with those of the inverse 
finite-size susceptibilities rff 1 . To this end we use the RG 
flow to map the original renormalizcd theory to one cor- 
responding to the choice I = (/x£oo) _1 of the scale param- 
eter. The running coupling constant u(£) can be replaced 
by the fixed-point value u* = 3e/(n + 8) + 0(e 2 ) at the 
expense of neglecting corrections to scaling ~ u(£) — u*. 
The running temperature variable f(l/^^ oc ) is exactly 
unity (at the required first order in u*, when r > 0). 

As straightforward consequences of Eqs. (4.31) and 
(4.32) we thus obtain 



R(pe r ) ( L ) = l + e^±| 167r2 ^ L2 ) +0 (e) 



and 

P>(sp,sp)^ 



1 + 6 



2 ^ 2 Q 4 , 2 (4L 2 )-7rL 3 



(4.46) 



n + 8 



L 4 



+ o(e). 
(4.47) 

Using the asymptotic forms (Dl) and (D5) of 0^2 (4L 2 ) 
for small and large L, one sees that these results are in 
conformity with the limiting behavior (4.36). The ampli- 
tudes Pq+ are found to be 



p^ r >=4p^+o(e) = e^ 



2 2^ 2 



8 3 



+ o(e). (4.48) 



The approach to the large-L limit R( p )(oo) = 1 is qual- 
itatively different for periodic and sp-sp boundary con- 
ditions: it is of an exponential and algebraic form in the 
first and latter cases, respectively. 

In Fig. 4 we have plotted the extrapolations to d = 3 
of the 0(e) results (4.46) and (4.47), obtained by set- 
ting e = 1, for the one-component case n = 1. It re- 
veals another important difference: The extrapolation 
R(p CT )(|_)| e=ri=1 remains positive for all L > 0, reassuring 
us thus that the theory is consistent in that the disor- 
dered state about which we expanded satisfies this nec- 
essary stability condition. By contrast, the extrapolation 
R(s P ,sp)( L )| e=n=1 b ecomes negative for 0.42 < L < 0.93. 
When extrapolated to d = 3 in this naive manner, the 
theory thus yields a violation of stability of the disordered 
state in this range of parameters. 

It is to be emphasized that this is a problem already 
for KD's original extrapolations of their e-expansion re- 
sults for the Casimir effect. As we shall see below, in our 
reformulated field theory it will show up in an even more 
exposed fashion. Note, however, that negative values of 
the 0(e) result for R( s p- s p)(L) are encountered only for 
values of e > 0.8265. This is illustrated in Fig. 5, where 
r( s p> s p)(L)|„ =1 i s plotted for several different values of e. 
It is conceivable, although not at all guaranteed, that ex- 
trapolations based on perturbative calculations to higher 
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FIG. 4: Plots of the scaling functions R (pcr) (L) and r( s p- s p)(L) 
for n = 1 and d = 3, as obtained from Eqs. (4.46) and (4.47) 
by setting e = n = 1. 
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FIG. 5: Plots of O(e) result for the scaling function R< s p- s p)(L) 
for n = 1 and the indicated values e = 0.1, 0.5, 0.8265, and 1. 



orders will yield positive definite functions R( sp,sp )(L). As 
already remarked above, we believe that in systematic 
studies of the stability of the disordered phase, besides 
temperature, the surface enhancement variables C\ and 
C2 should be allowed to vary — a difficult task, which is 
beyond the scope of our present analysis. 



E. Scaling functions of the residual free energies 

To determine the free-energy scaling functions (L), 
we start with the decompositions 



fies,R( L 'i T ' U > A 4 ) — //res./?. + /< 



Tp.rcs.R ' J <p : res,R 



(4.49) 



and 



eW(L) = e< p) (L) + e(f)(L). (4.50) 
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analogous to Eq. (4.2). We now substitute the pertur- 
bative expressions (4.18), (4.19), and (4.35) for 



and f^l^, insert Eqs. (4.31) and (4.32) for r K ^> and 
r (sp,sp)^ together w jth their common large-!/ limit 



(per) 



T + T ■ 



e/2 _ x 



6 W 2 e/2 



+ 0(u 2 ) (4.51) 



for the inverse bulk susceptibility r^. This yields the 

(truncated) series-expansion results for f^J R (L;T,u,fi) 
on which our subsequent analysis is based. We now com- 
bine them with the RG, proceeding along the lines ex- 
plained and followed above. 

The functions ©if^ have conventional expansions in 
integer powers of e, which to first order in e follow directly 
from Eqs. (4.18) and (4.19). Our results are 



0i PCr) (L) = -n 



4tt [Q 6 , 2 (L 2 ) -ei? 6 , 2 (L 2 )] 



L 2 



+ ne n + 8 
+ 0(e 2 ) 



2 [L 3 + 8^Q 4 , 2 (L 2 )]' 



8L 4 



(4.52) 



and 



q( s P. s p) 



(L) 



= — nn 



+ ne 



+ 0(e 2 
where i? 6-2 is defined by 

RdAr) 



(l + eln2)Q 6 , 2 (4L 2 )-ei? 6 , 2 (4L 2 ) 
16 L 2 

n + 2 ^Q 4 , 2 (4 L 2 ) [2L 3 + u Q 4 , 2 (4L 2 )] 
n + 8 

,2\ 



32 L 4 



dQdA r ) 
dd 



(4.53) 



(4.54) 



Inspection of KD's work reveals that the non-zero 
mode part of their 0(e) expression for 0( s p< s p) coincides 
with their result for @( D ' D \ By consistency, the latter 
should agree with our result (4.53) for 



(sp,sp) 



This is 

indeed the case, as can easily be verified by comparison, 
using the relation 



#6,2 (r) 



32^ 3 l\° E 3 
fl3/2,i (Vr/2) 



In"- 

7T 



53/2,0 (Vr/2) 



implied by Eq. (3.26). 

A consistency check can also be made for 0^ pcr ^ by 
noting that the contribution produced by the L 3 term in 
[. . .] 2 of Eq. (4.52) corresponds to the subtracted fco = 
part. Thus, by dropping it, we should recover KD's 
result for 0( per ) given in the third line of their equations 
(6.13). Confirming this is again straightforward by virtue 
of Eq. (4.55). 



We stress that unlike the full scaling functions O^ (L), 
their non-zero mode parts 0^ (L) do not in general decay 
exponentially as L — > oo and should not be expected to 
have this property. This is because the zero-mode pieces 
projected out involve contributions to the residual free 
energy density / ros ,i? that decay as 1/L. These imply 
contributions to ©^(L) that vary as |_ d ~ 2 in the large-L 
limit. Inspection of our result (4.52) shows that the 0(e) 
term of 0^°^ indeed grows as L 2 . By contrast, our 0(e) 
result (4.53) for 0^ p ' sp ) i s SC en to decay exponentially for 
large L because both the functions Qd.i and Rd,2 do so [cf. 
Eq. (D4)]. The absence of an analogous 0(e) contribution 
~ L 2 to 0^ p,sp ' i s due to the cancellation of the two 
terms of /^p] in Eq. (4.12) proportional to A d _ 1 /L and 
Bd/L, respectively. Of course, if such cancellation did 
not occur then the above-mentioned equality of 0^ ip ' sp '' 
with 0^f' £> '' to first order in e would be impossible. 

We next turn to the computation of the functions 0^ . 
This is a considerably more subtle problem, which re- 
quires care. It should be clear that we must not sim- 
ply expand in powers of e. The small- L behavior of 
the scaling functions O^ should be compatible with 
the behavior found for r = in Ref. [19] and hence 
yield the contributions ~ e 3 / 2 to the Casimir ampli- 
tudes. The mechanism by which this happens is that 
the inverse susceptibilities r^\r,u*) approach nonzero 
limits r* L (0,u*) = 0(u*) as r — > when L < oo. 
The 0(e 3 / 2 ) terms then result from the contributions 



„(d-l)/2 



3/2+0(e) 



On the other hand, if we expand in powers of e, tak- 
ing L (i.e., t) to be positive, then KD's series-expansion 
results to order e still ought to be recovered. 

Substitution of the respective one-loop results (4.31) 

and (4.32) for rf^ in the zero-mode free-energy contri- 
bution (4.35), in conjunction with Eq. (2.24), yields 



)] 3/2 } 



(4.56) 



where R( p )(L) represents the respective 0(e) expression 
for these scaling functions given in Eqs. (4.46) and (4.47). 

In the case of periodic boundary conditions, which 
we consider first, the combination of Eqs. (4.46), (4.50), 
(4.52), and (4.56) leads to 



(4.55) e( pcr )(L) 



= — n 



+ ne 



4^[0 6 , 2 (L 2 )-ei? 6 , 2 (L 2 )] 
L 2 

n + 2 [L 3 + 8^Q 4 , 2 (L 2 )] 2 -L 6 



n + 8 
nL 3 



127T 

+ o(e 3 / 2 ) 



1 + e 



8L 4 

n + 2 16tt 2 Q 4:2 (L 2 )^ 3/2 ' 



L 4 



(4.57) 
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This result has the following properties: 

(i) Upon expanding it to first order in e [ i.e., the term 
[. . .] 3 / 2 in Eq. (4.57)] when L = 0, one recovers KD's 
result. 

(ii) The limiting value e( por )(0) agrees with our 
(9(e 3 / 2 ) result for 



A (P er ) 



nir 2 nir 2 e r 
"90" + 180 



1 - C E - In 7T + 



2C'(4) 
C(4) 



5n + 2 



2n + 
+ 0(e 2 ) 



7T.7T 2 / n 



9V6 V 



n + 8 



3/2 



e 3/2 



(4.58) 



in Ref. [19]. 

(hi) The small-L behavior of G)( pcr )(L) differs from the 
requested one specified in Eq. (4.37) by terms oc e 3//2 L; 
we have 



0(p«)(L) 



A (pe 



+ 



e 3/2 L + 0(L 2 



/n + 2^ 3/2 
l^o ~ 2V6\n + 8 / 

(4.59) 

In KD's result the term linear in L that is at variance 
with the limiting form (4.37) is of first order in e; here it 



is of the 



rder e 3 / 2 to which we determined A 



(per) 



(iv) The large- L asymptotic behavior of 0( pcr )(L) is in 
conformity with Eq. (4.42), just as KD's result is. 

It is gratifying that our result has the properties (i), 
(ii), and (iv). On the other hand, it still does not fully 
comply with the small-L form (4.37) dictated by the ana- 
lyticity of the total finite-size free energy at T C:OC , though 
the violations now occur at the corresponding higher or- 
der e 3 / 2 . 

In Fig. 6 our result for the scaling function 6^ per ^(L) 
with n = 1 and d = 3, obtained by setting e = 1 in 
Eq. (4.57), is plotted and compared with its analog for 
KD's e-expansion result. The minimum in KD's extrap- 
olation result appears to be due to the inadequate han- 
dling of the zero-mode contributions. Our extrapolation 
gives a monotonic behavior at small L, which agrees bet- 
ter both with the Monte Carlo data of Ref. [23] as well 
as with improved, more recent ones [25, 26]. 

In Fig. 7 analogous extrapolations to d = 3 of the 
scaling functions for n = 2, n = 3, and n = oo are 
displayed, along with the exact spherical-model result for 
d = 3. The comparison with the extrapolations based on 
KD's 0(e) results displayed in Fig. 1 indicates, on the 
one hand, that the extrapolations for given n oscillate as 
the order of the series expansion is increased and, on the 
other hand, that the variations with order are the bigger 
the larger n is. 

Next, we consider the case of sp-sp boundary condi- 
tions. In discussing extrapolations to d = 3 dimensions, 
we shall restrict ourselves to the n = 1 component case. 
The reason should be clear: Only when n = 1 is a mul- 
ticritical point expected to occur at T CiOC and a finite 
enhancement of the surface interaction constants [45-47] . 




FIG. 6: Extrapolations to d — 3 of the scaling function 
9 (pcr) (L) for n = l, obtained by setting e = 1 in Eq. (4.57) 
and KD's original O(e) result, respectively. 




FIG. 7: Extrapolations to d — 3 of the scaling function 
Q(p er )(L) for n = 2, n = 3, and n = oo, obtained by set- 
ting e = 1 in Eq. (4.57). For comparison, the spherical-model 
result for d — 3 [5, 20], which is exact for n = oo, is also 
shown. 



A first problem was encountered in our investigation 
of the inverse finite-size susceptibility r £ ip ' sp ' ) : Our one- 
loop result for the scaling function R( s p> s p)(|_, n=l, e=l) 
becomes negative for 0.42 < L < 0.93. Clearly, con- 
vincing predictions for the scaling functions 0^ sp ' sp ^(L) 
in d = 3 dimensions must also fulfill necessary stability 
conditions such as the positive definitencss of r^. Thus 
the violation of this stability criterion of the disordered 
state is a problem even for extrapolations to d = 3 of 
KD's original 0(e) results. In our result given by the 
combination of Eqs. (4.47), (4.50), (4.53), and (4.56) it 
manifests itself in an obvious, striking manner: For val- 
ues of the scaling variable L in the mentioned interval, 
the extrapolation to d = 3 would yield complex numbers 
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for e( sp < sp )(L). 

A further problem occurs for large L: The contri- 
bution to e( sp ' sp )(L) originating from the term oc 1/L 
in curly brackets in Eq. (4.56) in conjunction with the 
part oc 1/L of R( sp ' sp ) produce a large-L behavior of the 
form 0(e 2 ) L + 0(e 3 ) L°. Thus, unless we subtract these 
asymptotic terms oc e 2 L and oc e 3 L°, our approximation 
for 0( s P< s P)(L) will not have a finite limit as L — > oo, 
and hence yield unacceptable results at d — 3 even 
in the regime L > 0.92 where the positivity condition 
R(s P , SP ) > o is satisfied. 

The combination of these two problems puts us in 
a bad position to suggest convincing extrapolations to 
d = 3. Let us, however, note some appealing properties 
the result given by Eqs. (4.47), (4.50), (4.53), and (4.56) 
has. All above properties (i)-(iii) of the small-L behavior 
hold just as in the case of periodic boundary conditions. 
That is, KD's O(e) results are recovered when the term 
[R(sp,s P )(|_)]3/2 in Eq ( 4 56 ) is expanded in powers of e. 

Second, the limiting value O( s p< s p)(0) reproduces the ex- 
pansion of the Casimir amplitude to order e 3 / 2 , 



. (sp.sp) 



+ 



nn 
1440 

5 n 



+ 



nir 2 e 
2880 



1-C E - ln(47r) + 



2C(4) 
C(4) 



2n + 8 
+ 0(e 2 ). 



mr 2 /n + 2 
72 V6 \n + 8 



3/2 



,3/2 



(4.60) 



Third, the term linear in L that violates the limiting 
form (4.37) is of order e 3 / 2 rather than linear in e. We 
have 

@(sp,sp) (L) M A (^P) + ^L(H±l) 3/2 e 3/2 L + 0{L 2y 

(4.61) 

Furthermore, the large-L behavior still is in accordance 
with Eq. (4.43) in the sense that the differences are of 
higher than first order in e. However, as already men- 
tioned, it would lead to extrapolations to d = 3 that 
grow ~ L in the limit L — + oo unless contributions of the 
form - 0(e 2 ) L + 0(e 3 ) L° are subtracted. 

In Fig. 8 we have plotted the extrapolated scaling 
function 6 (sp ' sp) (L) one obtains from Eqs. (4.47), (4.50), 
(4.53), and (4.56) upon setting e = 1, together with its 
analog (labeled KD) implied by the e-expansion result. 
The former function is depicted only for values L below 
the lower threshold ~ 0.42 beyond which the extrapo- 
lated scaling function R( s p< s p) of the inverse susceptibility 
becomes negative. We have refrained from displaying it 
(or appropriate modifications of it) for values larger than 
the upper positivity threshold ~ 0.93. In view of the 
0(e 2 ) corrections the result would require for large L to 
ensure its decay for L — > oo, we have no convincing rea- 
sons to expect such ad hoc modifications to yield much 
better results in this regime of L than the extrapolated e 
expansion. 

One might wonder whether the above problems could 
be avoided by a different choice of the free propagator 



G v in Eq. (4.22). For example, one might want to use 
one whose mass parameter is simply the sum of the free 
contribution f and the first-order perturbative correc- 
tion (3.20). We have in fact explored this possibility. 
It yields a modified scaling function (~>( s P' s p)(L) whose 
large-L behavior must be corrected by 0(e 2 ) contribu- 
tions to avoid unacceptable divergences. Once this is 
done, its extrapolation to e = 1 gives real values for all 
L. We refrain from displaying the results because we 
consider them unsatisfactory for two reasons. First of 
all, as explained before Eq. (4.22), we believe that the 
use of the inverse finite-size susceptibility r^ 1 as mass 
parameter is the more natural choice. Second, the fact 
that one is able to produce a well-defined extrapolated 
scaling function (~)( s P' s p)(L) does not cure the problem 
that the 0(e) result for the scaling function R( s p< s p)(L) of 
tl becomes negative when extrapolated to e = 1. Con- 
vincing improvments should yield meaningful extrapola- 
tion results for both the scaling function 0( s P' s p)(L) and 
tl within one and the same consistent approximation 
scheme. Evidently, further work is necessary to improve 
on the present unsatisfactory state of these results for 
sp-sp boundary conditions. 



On the other hand, the behavior of our results at small 
L may be expected to be superior to those based on the e- 
cxpansion. One indication is that, in the case of periodic 
boundary conditions, our results are in conformity with 
the exact solution in the large- n limit (see Sec. V). 
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FIG. 8: Extrapolations to d — 3 of the scaling function 
e (sp,s P ) f or n = i. The curve labelled "KD" corresponds to 
the O(e) results of Ref. [8] evaluated at e = 1; the other one 
is our result given by Eqs. (4.47), (4.50), (4.53), and (4.56), 
with e set to one. For further explanations, see main text. 
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V. COMPARISON WITH SPHERICAL-MODEL 
RESULTS FOR PERIODIC BOUNDARY 
CONDITIONS 

As is well known, for translation invariant systems re- 
sults that are exact in the limit n — > oo can be obtained 
from the exact solution of spherical models [48]. The 
self-consistent equations from which the scaling function 
0g^(L) for the spherical model with periodic bound- 
ary conditions must be determined can be found in the 
literature [5, 20-22]. Our aim here is to verify the 
consistency of our results for periodic boundary condi- 
tions with the exact solution of the spherical model for 
2 < d = A — e < 4. Making an analogous check for 
p = (sp, sp) is a much harder challenge and will not 
be attempted here. The reason is that the presence of 
surfaces in general destroys translation invariance per- 
pendicular to the boundary planes. The large-n limit 
of n-vector models on slabs with two parallel boundary 
planes 03 1 and 2$2 is known to correspond to a mod- 
ified spherical model involving separate constraints on 
the sums ^ - e laycr z S? of the squares of the spin vari- 
ables for each layer z [49]. The resulting self-consistent 
equations, while not difficult to determine, involve a z- 
dependent self-consistent pair interaction and so far have 
not been solved analytically. 

The exact solution for the spherical-model scaling func- 
tion e^° r) (L) = lim„^ oc e( por )(L)/n may be gleaned 
from Ref. [21], where this function was denoted as Y . 
It is given by 



47rQ d+2 , 2 (R ) 
Ro 



f(R d /2 -L d ) 



(5.1) 



where Rq = L 2 Rg^^C-) is a solution to 



2 Qd,2(Ro) 
R^ 



Ad (Rq 



(d-2)/2 _ , d-2 



L d - 2 ) . (5.2) 



The latter equation is easily solved for small e. Since 
Ad has a pole oc e _1 , the left-hand side starts to con- 
tribute at 0(e). One obtains 

Ro(L) = L 2 + e I^i^ +o(e) , (5.3) 
which becomes 



Ro(0) = e-7r 2 + o(e) 



(5.4) 



Eq. (5.1) can be combined as 

-A d d/2 47rQ d+ 2.2(Ro) 



R 



o 



-Ad 



d-1 



1 R (d-l)/2 
" K 



fe=0 



(-R ) fc , (5.5) 



where 



Ofc(d) = ^^2T- H fc + (! - rf )/ 2 ] C(l - d + 2k) . (5.6) 

Except for 02(d), which has a simple pole at d = 4, the 
coefficients afc(d) are regular at d = 4. We therefore sep- 
arate the contribution from the first term in the second 
line of Eq. (5.5) 

-A d -1 R (d-l)/2 

K o 



d-1 
-1 

127T 



L 2 + e 



16^ 2 Q 4 , 2 (L 2 ) 
L 2 



3/2 



+ ( £ 3/2 ),(5.7) 



where we substituted R by its expansion (5.3), and then 

expand the remaining contributions to ©sm^C-) m pow- 
ers of e. This gives 



e s P M r) (L) = 



12tt 



2 167T 2 Q4,2(L 2 ) 



L 2 + e- 



L 2 



3/2 



12tt 



+ e 



47r[Q 6 , 2 (L 2 )-ej? 6 , 2 (L 2 )] 
L 2 

8^ 2 Q 4 , 2 (L 2 ) 



L 2 



L Q 4 , 2 (L 2 ) 
4tt L 2 



+ o(e 3 / 2 ) . 



(5.8) 



The result agrees with the one for 0^ per ^(L)/n given in 
Eq. (4.57) if the factor (n + 2)/(n + 8) is replaced by its 
large-n limit (= 1). In particular, its value at L = 0, 



e 



(P cr )/ 



e 3 / 2 + 0(e 2 ) , 



7 „ , 2C'(4) 
2-^- ln7r + W 



9V6 



(5.9) 



coincides with the limit lim^oo cr ' /n of the expan- 
sion (4.58). The same holds for the coefficient of the term 
linear in L, for which we find 



rf©s P M r) (L) 



dL 



e 3/2 



L=0 



2V6 



+ ( £ 3 / 2 ) 



(5.10) 



which is consistent with Eq. (4.59). 



at the bulk critical point. Comparison of these results 
with ours for L 2 Rg I j^ r ' ) (L) contained in Eqs. (4.46) and 
(4.48) shows that the latter reduce to them in the limit 
n — > 00. 

Turning to OgM^C-)' we n °te that according to the 
representation (C9) of Qd+2,2, two contributions in 



VI. SUMMARY AND CONCLUDING 
REMARKS 

In this paper we have reconsidered the use of renor- 
malized field theory near the upper critical bulk dimen- 
sion d* — 4 to the study of finite-size scaling in slabs of 
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finite thickness and the thermodynamic Casimir effect. 
In previous work [19] it had become clear that in those 
cases where the boundary conditions involve zero modes 
in Landau theory at the bulk critical point, the conven- 
tional RG-improved perturbation theory based on the e 
expansion becomes ill-defined at T CiOC due to infrared sin- 
gularities. This could be remedied by means of a reorga- 
nization of field theory, which revealed that noninteger 
powers such as e 3 / 2 appear in the small-e expansion. 

Our main aim here was to examine how the calculation 
of scaling functions describing the large length-scale be- 
havior of the residual free energy and the Casimir force 
near the bulk critical point can be reconciled with these 
findings, so that the results of Ref. [19] for T = T Cj00 are 
recovered in the appropriate limit. 

We were able to show that consistent scaling functions 
can indeed be obtained both for the case of periodic and 
sp-sp boundary conditions. It became clear that the ill— 
dcfinedness of the conventional e-expansion theory due 
to zero modes manifests itself already at two-loop order 
inasmuch as contributions found at this order were found 
to have no power-series expansion in e at T Ci00 since they 
vary <~ e 3 / 2 . 

In calculations of crossover scaling functions by means 
of RG-improved perturbation theory near an upper crit- 
ical dimension one usually is faced with the following 
problem. The RG commonly achieves the proper expo- 
nentiation of the infrared singularities only at the unsta- 
ble fixed point. However, it does not normally do this 
- at least, not automatically — for the modified sin- 
gularities that occur as the scaled crossover variable be- 
comes large. Knowledge about the corresponding asymp- 
totic behavior frequently is obtained from other sources, 
such as RG analyses of a different model or fixed point, 
or short-distance expansion. Representative examples 
are the calculation of the two-point correlation func- 
tion [50], the crossover at a bicritical point [51, 52], and 
the crossover from critical to Goldstone-mode behavior 
in isotropic ferromagnets [53, 54]. To obtain and ver- 
ify the correct singularities of the behavior to which the 
crossover occurs by means of the e expansion, it must be 
supplemented by appropriate assumptions, or preferably 
knowledge, about the respective asymptotic forms. In 
some cases it has even been possible to design RG pro- 
cedures that yield the correct asymptotic behaviors at 
both the unstable fixed point as well as the stable one to 
which the crossover occurs [52, 54], albeit with somewhat 
limited range of applicability and success. 

Similar problems evidently had to be expected in the 
study of the problems considered here — finite-size ef- 
fects and thermodynamic Casimir forces. However, the 
challenges are actually greater and the difficulties more 
severe. Ideally, one would like to have a theory that 
has the power to correctly treat the infrared singularities 
at both the bulk critical point as well as the film crit- 
ical point and moreover is capable of handling the cor- 
responding dimensional crossover. For reasons discussed 
at the end of Sec. Ill, such ambitious goals would be un- 



realistic for a theory based on an expansion about the 
upper critical dimension. We therefore set out to reach 
more modest goals, namely: to modify and correct the 
previous theory by an appropriate treatment of the zero 
mode in such a way that (i) RG-improved perturbation 
theory becomes well-defined for temperatures T > T Ct<x> , 

(ii) reasonable scaling functions result whose limiting be- 
havior complies with the theory's predictions directly at 
T c oo and can be extrapolated to d — 3 dimensions, and 

(iii) hence bring it into a state comparable to the one it 
has for the non-zero-mode boundary conditions p = ap, 
(D, D), and (D,sp). 

We feel that, on the whole, our results are encourag- 
ing, in particular, for the case of periodic boundary con- 
ditions, where besides achieving (i)-(iii), we were able to 
demonstrate consistency with the exact large-n solution. 
Moreover, the scaling function obtained by extrapolation 
to d — 3, at least in the one-component case, appears to 
agree reasonably well with Monte Carlo results [25, 26]. 

The case of sp-sp boundary conditions turned out to 
be more delicate. First of all, we found that the one-loop 
expression for the scaling function R( s P' s p)(L) of the in- 
verse finite-size susceptibility becomes negative in a small 
regime of L = L/£oo when e exceeds the value ~ 0.8265 
(see Figs. 4 and 5). This tells us that all extrapolations of 
free-energy scaling functions and Casimir forces to d = 3 
based on approximations which yield the same one-loop 
scaling function R( s p< s p)(L) are questionable, at least in 
the regime where the positivity condition R( s p> s p)(L) > 
is violated. This applies both to KD's original extrapo- 
lation and ours (see Fig. 8). 

Our investigation of this case also revealed another 
problem: Perturbative RG calculations do not necessar- 
ily yield the correct asymptotic large-L behavior, at least 
not automatically. This applies even for the conventional 
e expansion in cases where no zero mode is present inas- 
much as the algebraic prefactors <~ L^- 1 )/ 2 appearing 
in the asymptotic exponential behaviors of 0^ per ^(L) and 
6 (sp ' sp) (L) given in Eqs. (4.42) and (4.43), respectively, 
are obtained only in e-expanded form. However, it is 
more troublesome in the cases studied here, especially, 
for p = (sp,sp). The reason may be understood as fol- 
lows. On the one hand, we encountered powers of inverse 
finite-size susceptibilities we had to retain to ensure con- 
sistency with the behavior at T Ci(X) . On the other hand, 
by expanding other contributions in e, L-dependent terms 
of order e 2 and higher are dropped which may be needed 
to cancel similar L-dependent contributions originating 
from the unexpanded powers of R in order to avoid in- 
correct or even divergent large-L behavior of the scaling 
functions. 

A qualitative difference between periodic and sp-sp 
boundary conditions is that the latter involve, even in 
the semi-infinite case L = oo, both d- and (d — 1)- 
dimensional critical behavior, rather than just a dimen- 
sional crossover. With hindsight it is therefore perhaps 
not too surprising that the latter turned out to be the 
more difficult case. 
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As remarked earlier, special surface transitions are ex- 
pected to occur in three bulk dimensions only in the 
n = 1 case. When n > 1, anisotropic special transitions 
should be possible if the continuous 0(n) symmetry is 
broken by an appropriate easy-axis spin anisotropy at 
the surface [45-47]. This is because surface phases with 
long-range order should not be thcrmodynamically stable 
at temperatures T > T c oo , by analogy with the Mermin- 
Wagncr theorem [55]. However, the 0(2) case is excep- 
tional in that a surface phase with quasi-long-range or- 
der should be possible. In fact, recent Monte Carlo work 
[56, 57] indicated that the surface phase transition is of 
Kosterlitz-Thouless type. Thus a multicritical surface- 
bulk point at which the line of these surface transitions 
reaches T CiOC should exist as well [12], and was reported 
to be found in the cited Monte Carlo analyses. 

Since the lambda transition of Helium involves a (real- 
valued) two-component order parameter, this 0(2) case 
is of potential relevance for Casimir forces in confined 
liquid He. In the case of 3 He- 4 He mixtures in contact 
with a substrate (see, e.g., [58, 59]), 4 He usually gets en- 
riched near the wall and a superfluid surface film may 
form there. Since order-parameter correlations decay al- 
gebraically in it, the bulk transition in the presence of 
such a critical surface phase is reminiscent of the special 
transition. Whether the central issue we were concerned 
with in this work — the presence of zero modes in Landau 
theory — arises also in the study of the thermodynamic 
Casimir effect in such systems and what its consequences 
are remains to be seen. A proper analysis of this ques- 
tion requires generalizations of our model. To describe 
mixtures, a second density besides the order parameter 
is needed. In addition, care must be taken to ensure a 
proper description of the Kosterlitz-Thouless-like surface 
transition. 

The present work suggests extensions and complemen- 
tary work along several lines. The situation in the case 
of sp-sp boundary conditions is rather unsatisfactory. To 
improve it, it would be desirable to extend our analy- 
sis by allowing the surface enhancement variables Cj to 
vary. In fact, in order to clarify the effects of finite size 
on the phase diagram, and in turn resolve the issue in 
which range of parameters the disordered state is ther- 
modynamically stable, such a generalization appears to 
be unavoidable. An appealing other aspect of it would 
be that by varying the Cj , one could smoothly interpolate 
between the boundary conditions p — (D,D), (D, sp), 
and (sp, sp). 

In view of the great technical and conceptual difficul- 
ties one is faced with in such analytical approaches, we 
believe that careful checks of their predictions by alter- 
native means such as Monte Carlo simulations are ab- 
solutely necessary. For a long time detailed studies of 
the thermodynamic Casimir effect by this method ex- 
isted only for the case of periodic boundary conditions 
[3, 23, 60]. However, recently new simulation strategies 
for investigating this effect in lattice spin systems with 
free boundary conditions have been developed [24-26]. 



As a result, systematic numerical studies of the ther- 
modynamic Casimir effect under all sorts of interesting 
boundary conditions have become possible. 

On the side of analytical theories, it would be interest- 
ing to explore whether the present approach can be com- 
bined with existing RG approaches at fixed dimension 
d for the study of bulk and surface critical phenomena 
[28-30, 61]. Another important challenge is to develop 
reliable analytical approaches by which the Casimir ef- 
fect can be investigated below the bulk and film critical 
temperatures. Recent investigations of the ordered phase 
based on Landau theory or RG-improved Landau theory 
[62, 63] certainly should not remain the final word since 
they fail to give correct descriptions of the critical be- 
havior at both the bulk critical point as well as at even- 
tual film critical points. In addition, they are known to 
be sometimes even qualitatively wrong inasmuch as they 
may predict phases with long-range order that can be 
shown to be destroyed by thermal fluctuations. 
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APPENDIX A: COMPUTATION OF REQUIRED 
INTEGRALS 

The one- and two-loop Feynman integrals for in- 
volve the free ^-propagator (3.13) at coincident points 
x = x'. Substitution of the cigenfunctions (3.3) and 
(3.4) into Eq. (3.13) yields 



2 r(d—l) ^ | 



and 



(Al) 



,(sp,sp), _ 2 f (d 1} ^ cos 2 (^mz/L) 



(A2) 

In order to compute the integrals l[ p ^ (L; f) introduced 
in Eq. (3.19), we add and subtract a summand with m = 
and then use Poisson's summation formula (see, e.g., 
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Eq. (4.8.28) in [64]) 



with 



E /(-)= E / * 



f(t)e 



2nijt/a 



(A3) 



Recalling the definition (3.24) of the functions Qd,a and 
performing the required momentum integrals, one arrives 
at the results for i[ pct \L;t) and /j sP:Sp '(L; f) given in 
Eqs. (3.21) and (3.22). 

Turning to the calculation of I 2 p \L;t), we note that 
for p = per we have 



lt l) (L;f)=[ir i, (L;r)} 



(P cr ) / 



(A4) 



as a consequence of translation invariance along the z di- 
rection. 

To compute l2 P ' sp \ we use the representation 



r(sp,sp) 



dz 

T 



[G%r\x-,x\r) 



- Ad-1) Ad-l) ^ 
cos 2 (fc m z)cos 2 (fc m /z) 



n 2 



and the fact that 



4 f ^ cos 2 (fc m z) cos 2 (fc m /z) = 1 + i (5„ 
Jo L * 



(A6) 



for k m = mit/L and m, m' ^ 0. Upon inserting the latter 
result into Eq. (A5), one arrives at 



r (sp : sp) 



(L;f) = [l[ sp > sp \L;r 



1 OO 

+ ^4 1 Eft+^" 3 . ( A7 ) 



m— 1 



where 4 was defined in Eq. (3.23). 

In order to evaluate the second term on the right-hand 
side of this equation, we employ the analytical continu- 
ation of the Epstein-Hurwitz zeta function discussed in 
appendix A of Ref. [65], namely 



Ec? 2 +« 2 r 

.?=i 

i 



T(s - 1/2) 



2a 2s 2a 2s - 1 r(s) 

OO 

^m s - 1 / 2 X (1 _ 2s)/2 (2™) 



(iraY 1 - 2 *)/ 2 
, s^l/2. (A8) 



A straightforward calculation then yields 



r(sp,sp) 



(L;f) = [l[ sp ' sp \L;r)Y 



A 



2d-4 orf_3 . °d-5/2 



+ 



2L 

Q 2 d-3, 2 (4f£ 2 ) 

22d-4 j-£2d-3 



(A9) 



8T(3 - d) T 2 [{d - l)/2] cos 2 (d7r/2) 



(A10) 



Aside from /j p '(i;r) (j = 1,2), we also need to cal- 
culate the integrals J^ p \L;t) introduced in Eq. (4.6) for 
p = per and p = (sp,sp). To this end we insert our 

above results for l[ p ^ into Eq. (4.6) and use the property 
(4.17), obtaining 



j( por )(£;r) 



2 A 



Ld-l 



±±f(d-l)/2 



2A d 



.d/2 



UA+ 2 (A11) 



and 



j(BP,B P )( i; ^ = j(p«)(2L;t) . 



(A12) 



The r — > limit on the right-hand side can be eval- 
uated in a straightforward fashion with the aid of the 
representation 



O (r) r(d+2) " V K ^< 2 ^ 



(2tt)<V 2 ~ 



-2)/2 



(A13) 



for d — > d + 2, where K v {x) is a modified Bessel func- 
tion of the second kind, and their well-known asymptotic 
behavior 



K v (x) = 2 u - l T{v)x- v + 0{x z - v ) 

for v > (see, e.g., Eq. (8.446) of Ref. [66]). 
The result 



(A14) 



lim 87rQw(r) = 2^ r(d/2) ((d) (A15) 



r— >0 



can be expressed in terms of either one of the one-loop 
Casimir amplitude ^c\i] given in Eq. (4.7). Inserting it 
into Eqs. (All) and (A12) finally gives 



j(per) (L . f) _ 2A C,[1] 2A d _ 1 (d _ 1)/2 



(per) 



and 

j(«P.BP)(£,;f) 

respectively. 



24, d/2 87rQ d+2 , 2 (rL 2 ) 



2 a( s p. s p) , . 



n 



Ld-l 



2A d d/2 7rQ d+2 , 2 (4fL 2 ) 
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APPENDIX B: EVALUATION OF /^(L; 0) 



In this appendix we present the calculation of the 

r(p) I 



one-loop free-energy contributions f^ p ),JL;f) defined in 



Eq. (4.4) for f = 0. 

Upon applying Poisson's summation formula (A3), we 
obtain 

oo „(d) 

/£[!] {L;0) = f ( » C o } +nLj2 ™ S (kjL) In , 2 

j=l J Q=(v,k) 

(Bl) 

and 



(B2) 



where 



with 



and 



^V>0 ~ J''' + ■'i/>,s,0 



f » f (d \ 2 



(B3) 



(B4) 



/£,o ) =2/S;o P) =-f J^^P 2 - (B5) 

The expressions (Bl) and (B2) can be evaluated along 
lines similar to those followed in Sec. V. of Ref. [8]. This 
gives 



n£$Z / cos(fcjL)ln<7 2 = L-^-^A^}, (B6) 



and 



^ /"(d) 

nL^ / cos(2fcji)lng 2 = L-( <i - 1 )A^ p) , (B7) 
j=i Jq=(pM) 

where A^fj^ are the one-loop Casimir amplitudes of 
Eq. (4.7). ' 



APPENDIX C: SERIES REPRESENTATIONS OF 
THE FUNCTIONS Q d , a {r) 

In this appendix we wish to derive representations of 
the functions QdA r ) as generalized power series and to 
establish their relation (3.26) with the functions g a ,b(z) 
utilized by KD. 

Let us define the integral 



Id,a(k,y) = K d _i / dpp 
Jo 



n d-2 (fc 2 +P 2 ) (CT - 2 ^ 

y + k 2 + p 2 



(CI) 



where K& denotes the usual factor 

Ad) nl-d -d/2 

K - s l 'fi'i-^w (C2) 

Then the right-hand side of Eq. (3.24) can be written as 
QdAr) = ^|[^(0,l)-/ d+1 , ff (0,l)VF]r( d+CT - 5 )/ 2 
+ J2ldAk,r)\, (C3) 

k=£0 > 

where here and below the summation extends over 

all nonzero k € 2-/tZ. 

When k = 0, the evaluation of the integral Id,a(k, r) is 
straightforward, giving 



2 l-d n (3-d)/2 r (d+o--5)/2 

T[(d- l)/2] cos[7r(d + cr)/2] ' 



(C4) 



In the case of nonzero values of k, we Taylor expand 
in r about r = to obtain 



IdAk,r)=J2 I df(k,0) r -, 



(C5) 



where the required partial derivatives I^\k,r) 



d^Id,<r{k,r)/dr^ at r — are given by 

T (0, j)n „ n , _ (-iyr(j + 1) T\j + (5 - d- a) /2] 



( 47r )(d-l)/2 | & |2j+5-«i-a p(j + 2 - (7/2) ■ 

(C6) 

We now substitute the Taylor series (C5) into Eq. (C3) 
and interchange the summations over j and k in the last 
term. Recalling the series expansion 



OO 



(C7) 



71=1 



of Riemann's zeta function, one can perform the k- 
summations analytically to obtain 

£ iffik, 0) = 2 4°f(2v, 0) ((2j +5-d-a). (C8) 

fe#0 

Using this result together with Eq. (C4), one arrives at 
the representation 



QdAr) 



2 -d n (3-d)/2 



T[(d- l)/2] cos[7r(d + ff)/2 

2 -d-l Tj.l-d/2 



„(d+<r-3)/2 



+ 



(d+<r-2)/2 



r(rf/2)sin[7r((i + cr)/2] 

ff(d -i)/2 v [ r[7 + (3-d-^)/2] 
r(j + l-a/2) 

C(2j + 3-rf-a) 



( 27r )2j-(T+2 



(C9) 
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For general values of d and a, the first two terms on the 
right-hand side of Eq. (C9) have branch-cut singularities. 
Cauchy's ratio test shows that the remaining power series 
(3rd term) is absolutely convergent for complex y inside 
a circle of radius (2-7r) 2 . 

As is known from Ref. [21], functions Qd,a with non- 
integer values of a are encountered in the study of finite- 
size effects of systems with long-range interactions. From 
the series expansion (C9) the asymptotic behavior of the 
functions Qd,a(r) as r — > can be read off easily even 
for such general values of a. This representation may, 
of course, also be employed to compute the functions 
Qd,a{r) by numerical means for values of r inside the 
radius of convergence of the series. 

To establish the relation (3.26) between Qd,a and the 
functions g a ,o(z) [cf. Eq. (3.25)] employed by KD, it is 
convenient to use the expansion (A13). Substituting the 
Bessel functions in it by their integral representation 



K v {z) 



(£/2)"V*F 

r( 



-zt 



(CIO) 



The function Rd,2(j) has an expansion in modified 
Bessel functions analogous to Eq. (A13), which follows 
from it by differentiation with respect to d. Using it we 
have determined Re,2(r) by numerical evaluation. The 
result is depicted in Fig. 10. 



0.002 
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and interchanging the integration with the summation 
over j immediately gives Eq. (3.26). 



APPENDIX D: NUMERICAL RESULTS FOR 
AND PROPERTIES OF THE REQUIRED 
FUNCTIONS Q d ,<r(r) 

In the present work only functions Qd,a( r ) with the 
special value a = 2 are needed. The purpose of the 
present appendix is to present numerical results for these 
functions. 

The expansion (A13) of these functions in terms of 
modified Bessel functions lends itself well to numerical 
evaluation. Figure 9 shows plots of the functions Qa,2 (0 
and Qe,2 (r), which were numerically determined via this 
representation. 



0.032 




FIG. 9: Plots of the functions Qi,2{r) and Qe,2(r), obtained 
by numerical evaluation of the series expansion (A13). 



FIG. 10: Plot of the function Rs^{r) as numerically obtained 
from the series expansion that results from Eq. (A13) by dif- 
ferentiation with respect to d. 

The asymptotic small-r forms of these functions can 
be determined in a straightforward fashion from the rep- 
resentation (C9). One obtains 

Qi,2{r) _ 1 y/r r\nr 
r r~o 24 ~~ 8tt ~~ 32tt 2 

l-2C g + 21n(47r) r | C(3) „ 2 



32tt 2 

C(5) T 3 + 0(r 4 ), 



256tt 4 



20487T 6 



(Dl) 



Q 6)2 (r) 



7T r r 3 / 2 r 2 In r 
r^o 360 ~ 96^ + 48tt 2 + 256tt 3 

_ 3/2-^,n(4») r2 + 0(r , MD2) 



and 



r^0 



C E + ln(47r) - 240C'(-3) - 8/3 
720 

_ C7 g -2 + 24C , (-l) + ln(47r) 
192tt 7 
ln(r/7r) + C E - 8/3 
96tt 2 

+ 0(r 2 ln 2 r). 



.3/2 



(D3) 



Their asymptotic forms for large values of r follow from 
(cf. equation (B33) of Ref. [21]) 



Qd,2(r) 



„(d+l)/4 



r^oo 2(2ir)( d - 1 )/ 2 



l + 0(r-^ 2 )} , (D4) 
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jiving 

Q4, 2 (r) 



and 



r l/4 



r r^oo 2 (2tt) 3 / 2 



Qe,2(r) r 3 / 4 



1 + Ofr- 1 / 2 ) 



(D5) 



^6, 2 (r) 



r 3/4 



- , .. e -V^l n J_[l + ( r -l/2) 
r r^oo 8 (2-7r) 5 /2 4^2 L v 



r r^oo 2 (2tt) 5 /2 



l + 0(r- 1/2 )] , (D6) respectively. 
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